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We present GW based self-energy calculations for the state of coexisting spin-density wave and c?-wave su- 
perconductivity in a series of cuprate superconductors. In these systems, the spin resonance spectrum exhibits 
the typical 'hour-glass' form, whose upward and downward dispersion branches come from the gapped spin- 
wave and magnetic scattering of Cooper pairs, respectively We show that the crossover between these two 
different dispersion features leads to an abrupt change in slope in the quasiparticle self-energy, and hence the 
low-energy kink commences in the single-particle quasiparticle spectrum. The calculated electron-bosonic cou- 
pling strength agrees well with experimental data as a function of temperature, doping and material. The results 
demonstrate that the electronic correlations dominate the quasiparticle spectra of cuprates near the low-energy 
kink, suggesting a relatively smaller role for phonons in this energy range. 
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I. INTRODUCTION 

Fundamental information for quantifying the strength of 
correlation effects and their role in magnetic and supercon- 
ducting (SC) properties comes from the determination of the 
magnitude of quasiparticle renormalization, and its intrigu- 
ing relationship with the quasiparticle line- shape. In con- 
ventional superconductors, the most convincing demonstra- 
tion that electron-phonon coupling is responsible for both dis- 
persion anomalies and superconductivity was obtained from 
the quantitative correspondence between features in the elec- 
tronic tunnelling conductance and the phonon spectrum mea- 
sured by inelastic neutron scattering.^ If unconventional su- 
perconductivity is mediated by exotic bosons a similar cor- 
respondence should hold, even though details of the coupling 
between these excitations and the quasiparticle state may vary. 
In cuprates, pnictide and heavy-fermion superconductors, the 
emergence of the SC state near the end-point of static antifer- 
romagnetic (AFM) Neel order makes these materials among 
the most promising contenders to realize spin- fluctuation me- 
diated pairing.^ ^ In heavy fermion systems, coupling of the 
spin fluctuations to the quasi-localized rare-earth /-electrons 
washes out signatures of the dispersion anomaly and compli- 
cates the analysis. Similarly, in iron-based superconductors, 
the presence of multiple bands, multiple gaps and possibly 
multiple 'hot-spots' ^ makes such comparative study difficult 
and ambiguous. 

On the other hand, single band cuprates should provide 
a clean system for understanding the microscopic origin of 
characteristic quasiparticle anomalies in their two-particle 
spectra. In particular, the experimental observation of sud- 
den changes or 'kinks' in the quasiparticle dispersion near 
50 — 70meV in cuprate superconductors^ has raised the 
hope that the bosonic excitations responsible for these kinks 
might also mediate electron pairing in these materials. How- 
ever, the kink appears at an energy scale where phonons, 
polarons,^^ plasmons,^^ and spin-fluctuations'^ can all con- 
tribute to its origin. In fact, significant isotope effects,^^ lat- 
tice coupling, and charge and spin fluctuations^^ ^ have all 



been reported in the kink energy range. Furthermore, poorly 
understood phenomena associated with underdoping, such as 
the normal state pseudogap, competition of superconductivity 
with incommensurate stripe physics, and nematic order which 
opens a gap over this energy scale make this problem more 
complex. The fundamental problem underlying this complex- 
ity comes from a lack of consistency between various spectro- 
scopies which are sensitive to different aspects of the ground 
state, which has hitherto precluded the development of a com- 
monly accepted theory. 

Here we go beyond these earlier postulates by perform- 
ing realistic calculations of the electronic excitation spectra 
including all channels of spin and charge degrees of free- 
dom in the coexisting spin-density wave (SDW) and su- 
perconducting (SC) ground state. 21 The calculations repro- 
duce the typical hourglass shape of the magnetic suscepti- 
bility, as revealed by inelastic neutron scattering. We find 
a new interpretation for this form, with the lower branch 
associated with scattering of the Bogolyubov quasiparticles 
and the upper branch corresponding to a gapped spin wave 
spectrum. These two branches meet to form a magnetic 
resonance peak at a material specific energy. In turn this 
peak interacts with the electronic dispersion via the GW self 
energy, causing an abrupt break in the dispersion known 
as the 'low-energy kink' (LEK). These low-energy bosons 
are qualitatively distinct from the strongly-correlated param- 
agnons responsible for the 'high-energy kink'.^^ We present 
the doping, temperature, and momentum dependences of the 
LEK for BisSrsCaCusOg (Bi2212), La2-:cSr^Cu04 (LSCO) 
and Nd2-a:Cea^Cu04 (NCCO). Significantly, the present re- 
sults are consistent with the bosonic spectra found earlier 
in neutron scattering, resonant inelastic x-ray spectroscopy 
(RIXS), Raman scattering, and studies of the optical 'glue' 
function, as well as with the excitations responsible for the 
high-energy kink seen in angle-resolved photoemission spec- 
troscopy (ARPES). 

This paper is organized as follows. In Sec. II, we describe 
the calculation of the GW self-energy due to spin and charge 
fluctuations in the state of coexisting SDW and d-wave super- 
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conductivity. In Sec. IIIA, we describe the microscopic ori- 
gin of the LEK. The calculated single-particle spectra in the 
LEK region are compared with experimental data in Sec. IIIB. 
The momentum, temperature and doping dependence of the 
LEK are given in Sec. IIIC. The corresponding values of the 
electron-boson coupling constant are compared with exper- 
iments in Sec. HID. Finally, we conclude in Sec. IV. Some 
technical details are found in the Appendices. 



II. FORMALISM 



Our starting Hamiltonian with competing Hubbard interac- 
tion and d—SC order is^^ 



k,k' 



(1) 



k,k' 



where cj^ ^(ck,cr) is the electronic creation (destruction) oper- 
ator with momentum k and spin a = ±, ek is the free particle 
dispersion, taken from a tight-binding parametrization of the 
first-principles band structure with no adjustable parameters, 
[obtained values are listed in Table |lj, and e^? is the chemical 
potential. The quadratic terms are expanded within Hartree- 
Fock formalism, and the wave SC gap is calculated using 

BCS formalism as = Vgk ^^Zk' 9k' {c]^' |C'_k',|)' where 
the d-wave structure factor is gk = cos kx(i — cos kya. The 
average is taken over the ground state with combined SC and 
spin-density-wave (SDW) order. Here the pairing interaction 
V is taken to be a momentum independent parameter which 
gives the experimental value of the SC gap A at T = and 
V(k, k') = VgkQk'- Similarly, the pseudogap is taken as 
U S, where S is the self-consistent mean-field SDW order pa- 
rameter S = (^<^k+Q cr^k,cr) at the SDW nesting vec- 
tor Q = (7r,7r). With this, the Hamiltonian in Eq. [T]can be 
diagonalized straightforwardly and the resulting quasiparticle 
dispersion consists of upper (u = +) and lower (v = —) mag- 
netic bands (U/LMB) further split by superconductivity: 



^k = ±((^r)' 
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Here 



+ vEok is the quasiparticle dispersion in the 
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non-superconducting SDW state, £^ok = [ (^k ) ^)^] 
and^± = [a±a+Q]/2. 

The unit cell doubling in the SDW state causes the correla- 
tion functions, such as the Lindhard susceptibility, to become 
tensors, with off diagonal terms in momentum space repre- 
sentation associated with umklapp processes at Q^^ In the 
SDW state, charge- and longitudinal susceptibilities become 
coupled at finite doping^^ . In common practice the transverse, 
longitudinal spin and charge susceptibilities are denoted as 



7 X^^ and respectively. We collect all the terms into a 
single notation as where a = a gives the charge and lon- 
gitudinal components and a = —a stands for the transverse 
component. The noninteracting Lindhard susceptibility in the 
SDW-BCS case is a 4 x 4 matrix whose components are^ 

Xr/(q,^) = Gi5(k,cr,icJn)Gsj(k + q,a,icJn 



(3) 



^ E ^'',^J E cZaZ'{K^.uj). (4) 



We obtain Eq. |4] from Eq. |3] after performing the Matsubara 
summation over n. G is the 4x4 single-particle Green's 
function constructed from Eq. [T] in the Nambu space. The 
summation indices viv') = ± gives upper and lowe magnetic 
bands, respectively. Here, the coherence factor due to SDW 
order in the particle-hole channel is 

^ \ ^Ok^Ok+q J 
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The SC coherence factors are 
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Lastly the index m represents the summation over three polar- 
ization bubbles related to the quasiparticle scattering (m = 1), 
quasiparticle pair creation (m = 2) and pair annihilation 
(m = 3), as defined by 



X^;^/(k,q,cj) = T 



+ + 



u ■ 



(7) 
(8) 



It is interesting to observe that is the particle-hole scat- 
tering term which becomes gapped in the SC state while the 
pair scattering terms, x^'^, only contribute in this region. We 
show below that this crossover from the normal to the super- 
conducting state gives rise to the low-energy kink. 

In the normal state, the 2 x 2 RPA susceptibility is obtained 
from the standard formula^^ 
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, [1 + aaUxMi^, ^)] X?f (q, uj) + U [x?f (q, a;)] ' 

XrPA — 1 1 1 9 5 

[1 - f/x?f (q, oj)] [1 + <Ta[/xf2"(q, OJ)] + aa [UxHi^, uj)] 

. , . [1 - [/xff (q, ^)] xfl(q, ^) + ^ [xff (q, ^)] ' 

[1 - c/x?f (q, '^)] [1 + ^^c/xf2"(q, ^)] + [f^x?2"(q, w)] 
xfpA 12/21 (q.'^) = xEISki^^ („) 

' / [1 - f/xff (q, w)] [1 + <Ta[/xf|(q, w)] + [f/xff (q, c.)] 



In the longitudinal and charge channel (a = a, denoted by 
subscript 'cz'), the RPA corrections do not introduce any new 
poles and thus all the normal state structure lies above the 
charge gap in the particle-hole continuum. Along the trans- 
verse direction (a = —a, denoted by subscript 't'), a linear 
spin-wave dispersion develops in the normal state which ex- 
tends to zero energy at QP The necessary condition to yield 
a gapless Goldstone mode is that Eqs. 9pT reduce to the self- 
consistent SDW order parameter at g = Q, which is indeed 
the case in the normal state. 

In the SC state, this zero energy spin-wave shifts to = 
I Akj^ I + I Akj^+Q |, due to theparticle-particle (and hole-hole) 
scattering terms x^'^ in Eq.js] These terms have finite inten- 
sity only if the SC gap changes sign at the 'hot-spot' Qp3^ 
see Eq. [6] Above the SC gap, the spin-wave term coming 
from EqT]7] is turned on. We show below that this crossover 
produces the characteristic hour-glass phenomenology and is 
ultimately responsible for the LEK. 

In a GW formalism, the W describes the dressing of the 
Green's function G by electron-hole pairs acting as bosonic 
modes. In our purely electronic model we calculate the lat- 
ter, which couple to the electrons causing an increase in the 
low-energy electronic mass and a shortening of the lifetime 
at higher energies. These effects are described in terms of 
a complex self-energy S which is calculated within the GW 
approximation to Feynman-Dy son perturbation ihtorp^ as 

G(k + q, cr', zcJn, ^p)r(k, q, zcj^, ^p)Xi'(q, ^p), (12) 

Here T is the vertex correction, modeled using Ward's iden- 
tity. Since the bosons are 'built' from fermions, self- 
consistency is required in calculating the self-energy, as in the 
quasiparticle-Giy approximation developed in Refs. |2T] [23l 
and|29| 



III. RESULTS 

A. Microscopic origin of low-energy kink 

In Figures [TJa) and (b), we present the imaginary part of 
the transverse and longitudinal + charge channels of the RPA 
susceptibilities (denoted by x'U Xc^' respectively) near 
Q = (tt, tt) along the diagonal direction, for LSCO at one 
representative doping x = 0.16. The momentum integrated 
values are given in Fig. [TJc). More extensive data for the 
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FIG. 1. (Color online) (a) Transverse susceptibility Xt is plotted 
along the diagonal direction centered at Q = (tv^tv) for LSCO at 
X = 0.16. Symbols are corresponding neutron data for LSCO at 
X — 0.16 (circles)^^ and LBCO at x = 1/8 (squares)^^. (b) Same as 
(a) but in the longitudinal -i- charge channel for LSCO at x = 0.16. 
(c) x" at Q as a function of energy, (d) Corresponding quasipar- 
ticle self-energy for both transverve (Ht) and longitudinal -i- charge 
channel (T^cz)- The magenta and cyan shadings differentiate the SC 
scattering and spin- wave spectra. The adjacent inset expands the en- 
ergy scale to clarify the nature of crossover energy scales present in 
the corresponding left hand side figure. 



whole Brillouin zone and over a large energy range and also 
for Bi2212 and NCCO are presented in Fig. [6] below. As 
we have shown in our earlier calculations, the normal state 
spin-wave dispersion becomes massive near Q due to SC gap 
opening. Below this spin gap, the Bogolyubov scattering of 
the SC states on the Fermi surface gives an oppositely disper- 
sive branch, leading to an hour-glass dispersive feature. 

The real part of the computed self-energy is plotted of 
Fig. [TJc) in the low-energy region. Above this energy scale 
the self-energy becomes dominated by normal state param- 
agnon contributions which are important for the study of the 
high-energy kink,^"^ but not in the present case. The origin of 
the low-energy kink is different. At low energy, shows a 
break in slope [see inset figure], rather than a peak, which cor- 
responds to the neutron mode discussed above.^^ This break 
in slope in occurring near 50 — 70meV depending on the 
material under study, leads to the LEK. 
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FIG. 2. (Color online) (al) Theoretical and (a2) experimental 
single-particle dispersion in the energy range of the LEK for LSCO. 
The blue and green dashed lines give the MDC peak positions for the 
theoretical and experimental spectra, while the non-interacting LDA 
dispersion is plotted as a red dashed line. (a3) The real part of the 
self-energy, calculated as the difference between the MDC peaks and 
the LDA dispersion. Dashed green line gives our computed imagi- 
nary part of the self-energy. Similar results are shown for Bi2212 
[Ref.[9l in (bl,b2,b3) and for NCCO [Ref. 10| in (cl,c2,c3). 



B. Single particle dispersion and the LEK 
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FIG. 3. (Color online) (al-a3) and (bl-b3): Momentum dependence 
of the LEK in optimally doped Bi2212 at both T < Tc and T > Tc 
The ARPES data (green line) are taken from Ref. 1 1 at the same mo- 
mentum. (cl-c3) Doping dependence of the LEK in Bi2212. Cor- 
responding experimental data^ ^^^^ are plotted by comparing Tc with 
universal dome feature as given in Appendix A. (dl-d3) Same as in 
(cl-c3) but for LSCO where the corresponding experimental data are 
obtained from Ref. 14 



energy uj.] The real part of the self-energy is typically de- 
termined as the difference between the MDC peaks and the 
bare LDA dispersion (black dashed lines). We plot the cor- 
responding experimental and theoretical estimates of SMn 
the right column of Fig. [2j which also displays the theoreti- 
cal values of T^" . For single layer systems, the LEK is around 
70 meV for LSCO, but 50 meV for NCCO both in theory and 
experiment while for Bi2212 our theory finds a larger value of 
the kink energy around 100 meV whereas experimental data 
show a kink near 70 meV. 



The self-energy dressed single-electron Green's function is 
G^^ = — S, where the Green' function and self-energies 
are 4x4 matrices as defined above The single particle spec- 
trum is then computed as A{k^ uj) = —lmGd,ii{k^ co')/7r, and 
plotted along the nodal directions in Fig. [2] for the three mate- 
rials LSCO, Bi2212, and NCCO.^^ The corresponding experi- 
mental datsPSniare given in the adjacent middle column. The- 
ory and experiment demonstrate good agreement both in the 
shape of the dispersion and in the associated spectral weight. 
Notably, our calculation neglects possible modulations of the 
spectral weight due to matrix element effects, which can be 
important in ARPES^"^, STM^^, inelastic light scattering^^, 
and other spectroscopic j^^^-^^l For ease in comparison, the 
dashed lines in both figures show the dispersion calculated as 
the position of the peaks in the theoretical (blue) and exper- 
imental (red) momentum distribution curves (MDCs). [The 
MDC is defined as a plot of A{k^uj) vs momentum k at fixed 



C. Momentum, temperature and doping dependence of LEK 

In Fig. |3] we study the momentum (k), temperature (T), 
doping and material dependence of the LEK. At each T 
and doping, both the SDW and SC gaps are evaluated self- 
consistently assuming T independent values of U and pairing 
interaction V. The calculated single-particle spectra yield a 
very good description of the ARPES data over the entire Bril- 
louin zone. Some discrepancies are noticeable in the antinodal 
region where the pseudogap and SC gaps are the largest. Here, 
our theory slightly overestimates the experimental kink, which 
is also the case when the pseudogap increases with underdop- 
ing (see lower panel). Also, in Bi2212 we have neglected the 
bilayer splitting of the Cu02 bands, which is largest in the 
antinodal direction. Note that the agreement with experiment 
would not be significantly improved by including a phonon 
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FIG. 4. (Color online) (a) Comparison of experimental and theoret- 
ical electron-boson coupling constant A as a function of doping and 
material. The horizontal scale is normalized to the optimal doping 
for each material, (b) Same as (a) but as a function of k and T for 
Bi2212. It is important to distinguish our obtained values of A with 
most of earlier data [see for example Ref. Q. We calculate A with 
respect to the actual LDA dispersion whose slope is doping and mo- 
mentum dependent. In common practice, A is obtained assuming a 
bare linear dispersion which is often taken to be doping and momen- 
tum independent, (c)-(d) Difference of the estimation of A between 
experiment and theory, using the same symbols as in the correspond- 
ing top panels. 



contribution to the kinkP^ 



D. Coupling constant 

Fig. |4] shows the corresponding electron-boson coupling 
constant A, calculated from the gradient of which in the 
low-energy region becomes II' (^k) = — A^k- A in princi- 
ple varies with u but weakly in the low-energy region below 
70 meV and thus we report the values of A at this energy. As 
shown in Fig.|4] the experiments and theory agree well in the 
optimal doping region for all three materials (see figure cap- 
tion for details). The difference between experimental and 
theoretical As becomes more prominent as one moves away 
from the optimal region. When this difference is positive, it is 
suggestive of the presence of an additional bosonic coupling, 
including phonons. Note that this additional component is al- 
ways small compared to the main electronic component that 
we have calculated. 



modes responsible for high-energy kink (not a crossover).!^ 
The associated electron-boson coupling constant A depends 
strongly on the slope of the lower dispersion branch of the 
hourglass, and thus on the FS which leads to a doping, tem- 
perature and material dependence of the coupling constant. 
The susceptibility energy scales are observed directly in the 
inelastic neutron scattering and optical glue measurements, 
and should also be seen in RIXS and Raman spectra with im- 
proved resolution. We argue that such a crossover feature is 
a general phenomenon of any coexisting SDW and unconven- 
tional SC state and thus should also be present in iron-based 
and heavy fermion superconductors. 



ACKNOWLEDGMENTS 

This work is supported by the U.S.D.O.E grant DE-FG02- 
07ER46352, and benefited from the allocation of supercom- 
puter time at NERSC and Northeastern University's Advanced 
Scientific Computation Center (ASCC). The work at Los 
Alamos National Laboratory is funded by US DOE, BES and 
LDRD. 



Appendix A: Tight-binding parameters and the interaction 
strength. 

We obtain the tight-binding parameters by fitting to the ma- 
terial specific LDA dispersion^^ as listed in Table |l] The 
self-energy renormalizes the LDA band to obtain the experi- 
mental dispersion. The interaction parameters, both Hubbard 
U s and superconducting pairing potential Vs, are determined 
from fits to other experiments - specifically the values of the 
superconducting and pseudogaps seen in ARPES, optical and 
STM etc measurements - and the resulting values are listed in 
Table [nl^^ The values of U match well with the calculated 
values of the screened [/, as shown in Fig. |5] Notably, we 
find that over the doping range x > 0.05 the bare U is univer- 
sal - the same for all materials and dopings. Finally, we find 
that our computational method also provides a good model for 
ARPES and optical properties of the cuprates.^I^ 



Material 


t 








Z 


NCCO [Ref. 401 


0.42 


-0.1 


0.065 


0.0075 


0.4 


LSCO [Ref. 40 1 


0.4195 


-0.0375 


0.018 


0.034 


0.3 


Bi2212 [Ref. 35| 


0.44 


-0.1 


0.05 


0.0 


0.4 



TABLE I. Tight-binding parameters obtained by fitting to the LDA 
band-structure with self-consistent renormalization factor Z (right). 



IV. CONCLUSION 

In summary, we have shown that the low-energy kink 
arises mainly from a crossover from Cooper pair scattering to 
spin- wave dispersion which is different from the paramagnon 



Appendix B: Material dependence of the susceptibility and the 
self-energy 

Calculated spectra of x" ^re presented in the first two 
columns of the Fig. [6] as a function of excitation energy along 
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Material 


Doping (x) 


Apg (meV) 


U/t 


Asc (meV) 


Pairing Potential 


TcK 


Z 






(Exp./Theory) 


(Theory) 


(Exp./Theory) 


V (meV) (Theory) 


Exp. (Theory) 




LSCO 


0.06 


150 


2.35 


6 


-93 


18 (27) 


0.5 


LSCO 


0.12 


120 


2.27 


11 


-63 


30 (48) 


0.48 


LSCO 


0.16 


63 


2.25 


15 


-51 


40 (85) 


0.45 


LSCO 


0.18 


43 


2.25 


13 


-35 


37 (75) 


0.43 


LSCO 


0.22 





2.25 


8 


-28 


26 (48) 


0.4 


NCCO 


0.15 


170 (athotspot) [Ref. 44| 


4.1 


5.5 rRef.l45l 


-83 


24 (31) rRef.l45l 


0.4 


Bi2212 


0.10 


113 


2.46 


15 


-11 


55 (85) 


0.45 


Bi2212 


0.12 


95 


2.42 


17.5 


-75 


65 (95) 


0.43 


Bi2212 


0.16 


75 


2.36 


20 


-67 


91 (115) 


0.4 


Bi2212 


0.19 


50 


2.36 


17.5 


-58 


70 (90) 


0.39 


Bi2212 


0.22 


25 


2.36 


12.5 


-50 


55 (75) 


0.38 



TABLE IL The value of U/t is chosen to reproduce the experimental pseudogap {Apg) whose values are presented here along the hot-spot 
direction in the electron doped case and the antinodal direction for hole doped cuprates LSCO and Bi2212. Similarly, the parameter value of 
pairing potential V is taken to reproduce the superconducting gap {Asc), whose maximum value along the antinodal direction is presented 
here. Our mean-field calculations overestimate the values of Tc, presumably due to the neglect of phase fluctuation^.Experimental data for 
LSCO are taken from Ref. 47 and for Bi2212 from Ref. 48. 



— LSCO (ARPES,Optical,RIXS) 




doping (x) 

FIG. 5. (Color online) Doping dependent values of U. The solid line 
gives the self-consistent values of [/ as a function of doping found 
in the QP-GW model^^ Symbols give the values used in the present 
calculations for various materials. 



high- symmetry directions in momentum space for the trans- 
verse (column 1) and longitudinal spin plus charge (column 2) 
channels; the three rows represent different materials, LSCO, 
Bi2212 and NCCO. Superimposed on the calculated spectra, 
we show the single magnon RIXS results for undoped LCC^, 
neutron data of the same sample^^, and RIXS data of undoped 
SCOCand NCO.i9 20 30 49 

We plot the corresponding momentum-averaged Eliash- 
berg, or 'glue' function in the third column of the Fig. [6] 
These functions display several peaks, but the most relevant 
ones for the present purpose are the two which are present 
in the energy scales of 300-400meV and ~ 70meV for hole 
doping and ~ SOOmeV and 50meV for electron doping. The 
high energy peak mainly stems from the strong intensities near 
(tt, 0) and lies in the 'waterfall' or high-energy kink region, 
while the low-energy peak arises from the magnetic resonance 
scattering around (tt, tt). 

The insets in column 3 of Fig.|6]show expanded views of the 



low-energy region of the a^F{uj) functions and compare them 
with available data from optical absorption^ ^ and ARPES 
(black lines)^^ In ARPES experiments, the a^F spectra are 
usually derived from maximum entropy calculations which 
tend to produce spectra consisting of a series of peaks, in- 
cluding the one at ~ 50 — 70meV. A similar multi-peak struc- 
ture is seen in the optical-glue function of NCCO (inset j^^. 
While both experiments find the feature at 50 — 70meV, 
there are additional peaks in the low-energy region not repro- 
duced in our calculation which may represent contributions 
due to phonons or other bosons. For NCCO, our theory pro- 
duces a kink rather than a sharp peak at ~ 50meV, which is 
consistent with the weaker strength of the LEK as compared 
to LSCO and Bi2212, see main Fig. [6] 

The real part of the computed self-energies for both spin 
channels and the charge channel are plotted in the last col- 
umn of Fig. [6] for all three materials. Here, again we see 
two energy scales of different characteristics. All the spin 
and charge components of are linear in the low-energy 
region— coming from the linear dispersion of the fluctuation 
spectrum along F (tt, 0) and F (7r/2,7r/2), shown in 
the left two columns— giving a total dispersion renormaliza- 
tion of the order of 2-3, consistent with experiments. Tj' at- 
tains a peak around 200meV [400 meV for electron doped] 
which means that the resulting electronic dispersion also un- 
dergoes a change in its renormalization behavior, yielding the 
high-energy kink. 

The LEK in T.' occurs at 70meV for LSCO and Bi2212 
and at 50meV for NCCO, consistent with ARPES data (black 
lines). The ARPES -derived usually shows a more pro- 
nounced peak at the low-energy kink, rather than a break-in- 
slope. This is partially due to the assumed form of the bare 
dispersion, which is taken as a straight line from the Fermi 
level to the dressed band at a high energy usually chosen at 
-300 meV, rather than the correct LDA band. Note that due 
to the Kramers-Kronig relation, the imaginary part of the self- 
energy increases gradually with energy and does not show any 
characteristic features at these two kink energies but shows a 
maximum at the energy where changes sign. Therefore, 
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FIG. 6. (Color online) The imaginary part of the susceptibility is plotted along the high- symmetry lines of the Brillouin zone for transverse 
[first column] and longitudinal plus charge channels [second column] for two hole doped (LSCO and Bi2212) and one electron doped cuprate 
(NCCO) near their optimal dopings. The results are compared with RIXS data for insulating LCO [blue from Ref 18], underdoped LSCO 
[black and cyan symbols from Ref. 18], and undoped SCOC [grey symbols]. Magenta and green symbols are neutron data for undoped 
and optimally doped LSCC^^"^^ which are compared with the LBCO data at 1/8 doping (deep green symbols)^^. The corresponding bosonic 
spectral weight a'^F{(jj) (see text) is compared for transverse (red) and the mixed state of longitudinal spin and charge (green) channels in the 
third column. The corresponding insets expand the low-energy part of the calculated a^F(u) — J d^Qx'iQi ^) compare them with 
data extracted from optical (cyan dashed lines) and ARPES spectra (black dashed line). The optical data are taken for nearly optimally doped 
LSCO and for Bi2212 from Refs. 31 and 32 and electron doped PCCO from Ref. 50 ARPES-derived data are for LSCO at x = 0.03^^ and 
overdoped Bi2212^^. The corresponding quasiparticle self-energy is presented in the right column. Black dashed lines are ARPES self-energy 
data from which a^F was extracted. 
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